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ABSTRACT 


After  defining  geomagnetic  noise  as  it  applies  to  MAD, 
the  geomagnetic  indices  currently  used  by  the  fleet  to 
predict  MAD  geomagnetic  noise  are  reviewed  to  determine 
their  actual  applicability.  The  current  indices  are  deter¬ 
mined  to  be  insufficient,  methods  are  proposed  for  estab¬ 
lishing  a  new  NAD  index,  and  a  developmental  NAD  index 
system  was  tested.  Geomagnetic  fluctuations  in  the  .04  to 
2.0  Hz  freguency  band  were  recorded  at  Monterey,  California, 
and  used  for  a  preliminary  test  of  the  proposed  MAD  index. 
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I.  SBfiaASSHISS  asirss 


A.  HISTORY  OP  GEO MAG NET ICS 

The  beginnings  of  the  study  of  geomagnetism  lie  back  in 
prehistory  when  magnetic  attraction  between  iron  and  certain 
minerals  was  first  observed.  Exactly  when  this  phenomenon 
was  first  noticed  is  not  known,  but  the  properties  of  magne¬ 
tite,  then  called  lodestone,  appeared  in  Greek  literature 
around  600  B.  C.  (Brennan  and  Davis).  [Ref.  1] 

Chapman  [Ref.  2]  indicates  that  the  directional  property 
of  magnets  was  known  and  used  in  Europe  prior  to  1200  A.  D. 
and  possibly  in  China  before  then.  E.  N.  Parker  (Ref.  3] 
noted  "It  is  an  interesting  fact  that  the  ancient  walls  of 
Peking  were  lined  up  with  magnetic  north  rather  ■‘■han 
geographic  north,  a  difference  at  that  time  of  about  lo.  we 
may  presume  that  the  surveyor  found  it  easier  to  work  with 
his  compass  needle  by  day  than  to  sight  on  the  pole  star  by 
night.”  This  property  allows  the  use  of  the  Earth's  geomag¬ 
netic  field  for  navigational  purposes. 

By  the  mid-fifteenth  century  it  was  determined  in  Europe 
that  the  magnetic  compass  does  not  point  to  true  north.  The 
angle  between  true  north  and  the  direction  indicated  by  the 
compass  is  now  known  as  magnetic  declination  by  the  geophy¬ 
sicist  and  as  variation  by  the  navigator. 

The  magnetic  field  dip,  or  magnetic  inclination  is  the 
angle,  in  a  vertical  plane,  between  the  horizontal  and  the 
direction  of  the  earth's  magnetic  field  vector.  It  was 
observed  in  1544  by  an  instrument  maker  in  Nuremberg  named 
Hartmann,  and  again  by  Robert  Norman  in  London  in  1581. 

These  discoveries  or  observations  gave  rise  to  the  study 
of  geomagnetics  as  a  specialty  which  gained  its  cornerstone 
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with  William  Gilbert  in  160  0-  After  comparing  his  experi¬ 
mental  results  with  the  previous  work  of  others  sur  as 
Norman,  Gilbert,  in  his  book,  "Jit  54311411# n  concluded  that 
Jftl&fifi  i£Sl  Sfit  aliens  li£Clfil£iS  (the  earth  globe 
itself  is  a  great  magnet)"  [ Bef .  2].  It  is  this  concept, 
that  the  earth  is  itself  a  magnet,  that  is  the  basis  of  the 
science  of  geoaagnetics. 

Gilbert  felt  that  the  Earth* 3  magnetism  must  remain 
constant  except  for  geological  changes,  but  it  was  soon 
determined  that  this  was  not  the  case.  A  'secular  varia¬ 
tion’  of  the  Earth's  was  found  to  exist. 

Shorter  term  changes  in  the  geomagnetic  field  were 
observed  and  it  was  eventually  realized  that  geomagnetism  is 
dynamic.  In  1722  George  Graham  discovered  that  daily,  or 
diurnal,  variations  exist  [Bef.  2].  During  the  early 
Nineteenth  century  magnetic  observatories  began  to  be  estab¬ 
lished  to  record  the  changes  in  the  geomagnetic  field  in  a 
systematic  fashion  (Knecht}  [Bef.  4]. 

B.  EARTH'S  MAGNETIC  FIELD 

i.  sanslilasals  °£  Hi  Ssaiaaasi&s 

There  are  various  ways  of  breaking  down  the  consti¬ 
tuting  parts  of  the  geomagnetic  field.  One  way  is  to  divide 
the  field  in  terms  of  distance  from  the  center  of  the  earth. 
Doing  this  yields  these  three  parts:  internal,  crustal,  and 
external  ( AFGL)  [Bef.  5].  The  internal  field  originates  in 
the  core  region  and  is  the  more  stable  field,  containing 
only  extremely  low  frequency  temporal  variations.  The 
crustal  (or  anomalous)  field  arises  from  modifications  made 
on  the  internal  field  by  materials  and  structures  in  the 
Earth's  crust.  These  variations  are  not  spatially  constant 
and  give  rise  to  some  of  what  is  known  as  geological 
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variations.  The  external  field  is  the  most  dynamic  and 
arises  from  many  sources  including  the  interaction  between 
the  solar  wind  and  the  terrestrial  magnetic  field. 

Another  way  of  describing  the  components  of  the 
geomagnetic  field  is  by  time  variation.  This  division  is 
accomplished  by  considering  that  part  of  the  field  which 
varies  with  periodicities  greater  than  about  one  year  as  the 
Siaail  £Asl&  and  what  is  left  as  the  yajlaiAail  £ASli 
(Knecht)  .  £Bef.  4] 

The  steady  field  consists  of  the  above  named 
internal  field,  also  referred  to  as  the  main  £iai£.  Slow 
variations  in  the  main  field  with  periods  of  years  or  longer 
are  referred  to  as  secular  lajiaiigUS* 

2.  jtedsis  a£  ika  aiiu  XAslA 

Various  models  of  the  main  geomagnetic  field  mate 
use  of  a  geocentric  aipola.  Gauss,  in  1839,  demonstrated 
that,  as  a  fairly  good  first  approximation  to  the  geomag¬ 
netic  field,  the  field  of  a  uniformly  magnetized  sphere  (for 
points  outside  the  sphere)  is  equivalent  to  the  field  of  a 
magnetic  dipole  located  at  the  center  of  the  sphere 
(Jacobs).  [Ref.  6] 

The  simplest  of  the  present  approximations  of  the 
geomagnetic  field  is  that  of  a  short  bar  magnet  or  dipole 
located  at  the  center  of  the  earth  with  an  axis  inclined 
approximately  11.5°  from  the  Earth's  axis  of  rotation.  The 
sense  of  the  field  lines  is  from  south  to  north  (Figure 
1.1)  . 

The  axis  of  this  dipole  intersects  the  earth  at  the 
geomagnetic  north  pole,  78.  5°N,  291. 0°E  in  geographic  coor¬ 

dinates,  and  at  the  geomagnetic  south  pole,  78.50S,  110.0°E. 
The  moment  of  the  geomagnetic  dipole  is  8.1  x  10*2  amp-m*. 
It  is  these  poles  that  are  used  to  define  the  geomagnetic 


Figure  1. 1  Dipole  ippe&rance  of  Geomagnetic  Field 


coordinate  system  (Knecfat)  [Ref  -  4].  The  geomagnetic  coor¬ 
dinate  system  is  a  spherical  polar  system  similar  to  the 
geographic  coordinate  system  with  a  geomagnetic  equator 
defined  90  degrees  away  from  either  geomagnetic  pole  in 
latitude-  This  tilled  geocentric  dipole  model  describes  the 
main  geomagnetic  field  to  an  accuracy  of  about  1055- 

In  1940  Chapman  and  Bartels  defined  an  off-center 
dipole  in  the  earth's  interior,  called  the  eccentric  dipole. 
This  dipole  is  displaced  0.  0685  earth  radii  (436  km)  in 
magnitude  from  the  center  and  in  the  direction  of  the  point 
15.6°N,  150. 9°S  (geographic  coordinates)  (Vestine)  [Ref.  7], 


The  intersections  of  the  eccentric  dipole  axis  at  the 
earth's  surface  are  81 • 0 ® 84,7®M  and  75.0°S#  120. 

(Figure  1.2)  (Hayaes)  [Ref.  8].  This  approximation  is  accu¬ 
rate  to  within  a  few  percent. 


The  field  has  additionally  bsen  modeled  to  an  accu¬ 
racy  of  about  lid  by  determi nining  Saussian  coefficients  by  a 
least-squares  fit  of  experimental  measurements  of  the 
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geomagnetic  field.  These  coefficients  are  used  in  a  spher¬ 
ical  harmonic  series  repra sen ting  the  scalar  potential  of 
the  field.  This  accuracy  implies  that  the  internal  contri¬ 
bution  to  the  total  main  field  is  at  least  on  the  order  of 
9  9%. 

The  International  Seonagnetic  Reference  Field  (IGHF) 
yields  values  which  differ  by  only  parts  per  thousand  from 
measured  values. 

3.  £&U£££fi  2t  *£21££&iU£  E«£l& 

There  are  various  elements  that  contribute  to  the 
geomagnetic  field,  some  external  to  the  earth's  surface  and 
some  internal.  As  previously  meationned,  the  external 
contributions  make  up  only  a  small  fraction  of  the  steady 
field,  playing  a  more  important  role  in  the  variation  £1^14. 

These  external  sources  Include  current  systems  in 
the  earth's  upper  atmosphere  affected  by  solar  electromag¬ 
netic  radiation  and  gravitation,  solar  corpuscular  radiation, 
or  the  interaction  of  solar  plasma  with  the  main  field,  and 
the  effect  of  the  solar  interplanetary  field.  [Ref.  4] 

Various  magnetic  surveys  of  the  world,  including 
those  conducted  at  ground  level,  by  airborne  instruments, 
and  by  satellite, have  pointed  to  the  fact  that  the  largest 
source  of  the  earth's  magnetic  field  is  internal  to  it. 
While  there  exists  residual  permanent:  magnetism  in  the 
earth's  crust,  this  cannot  be  the  principal  internal  source 
of  the  geomagnetic  field  due  to  temperature  and  material 
properties  known  to  exist  in  the  earth's  interior  (Nagata 
and  Ozima)  .  [Ref.  9  ] 

Permanent  magnetism  is  generated  by  microscopic 
electric  currents,  since  a  changing  electric  field  will 
generate  a  magnetic  field.  Another  way  *■  o  generate  a 
magnetic  field  is  by  the  motion  of  electric  charges  in  a 
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macroscopic  current.  Convective  motion  of  the  electrically 
conducting  fluid  core  of  the  earth,  resulting  in  a  macro¬ 
scopic  current  system,  is  considered  to  be  the  principal 
source  of  the  main  field. 

The  most  promising  present  theory  of  the  generation 
of  the  geomagnetic  field  is  that  of  some  sort  of  a  self-ex¬ 
citing  dynamo  system.  This  means  that  the  motion  of  a 
conductor,  such  as  the  molten  iron  in  the  earth's  core,  in  a 
magnetic  field  produces  a  current  which  in  turn  induces  a 
magnetic  field  in  support  of  the  original  magnetic  field 
[Ref.  4,  7].  K  very  simple  model  of  such  a  dynamo  is  shown 
in  Figure  1.3  . 


O 


Figure  1.3  Simple  Disk  Dynamo 


The  original  excitation  or  seed  field  nay  be  due  to 
an  external  field  line,  perhaps  from  the  solar  interplane¬ 
tary  magnetic  field.  This  original,  poloidal,  field  line  is 
wound  up  due  to  the  differential  rotation  (rotation  not 
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constant  with  latitude)  in  the  molten  core.  The  wound  up 
line  becomes  an  intense  azimuthal  field  which  is  carried 
outward  by  the  upwelling  associated  with  convection  and 
twisted  by  the  action  of  the  Coriolis  force.  The  twisting 
generates  a  helical  toroidal  field  which ,  by  outward  diffu¬ 
sion,  generates  the  eitocnally  observed  quasidipole  geomag¬ 
netic  field. 

The  combination  and  interaction  of  two  (see  Figure 
1.4)  or  more  dish  dynamos  can  also  can  also  explain  the 
reversal  of  the  geomagnetic  field  [Ref.  7,  9]. 


r 


Figure  1. 4  Twin  Dish  Dynamo 


Regional  anomalies  that  are  nondipolar  (do  not 
cor  form  to  the  dipole  field)  possibly  arise  from  eddy  circu¬ 
lations  in  the  outer  core  [  Hefi.  4]. 
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4  •  Magnetosphere 


The  magnetosphere  can  be  defined  to  be  that  region 
(see  figure  1.5)  occupied  by  the  geomagnetic  field  above  the 
ionosphere,  a  region  where  the  field  atongly  influences  the 
dynamics  of  ionized  gas  and  charged  particles  (Kern) 
C  Be  f.  10]. 

If  the  space  surrounding  the  earth  were  a  perfect 
vacuum,  the  marth's  magnetosphere  and  magnetic  field  might 
be  more  or  less  symmetric  and  extend  outward  until  it  merged 
with,  and  its  strength  became  insignificant  compared  to  the 
solar  ar.d  other  planetary  magnetic  fields.  This  turns  out 
not  to  be  the  case. 

Instead  of  being  an  island  in  a  perfect  vacuum,  the 
earth  encounters  a  continuous  flow  of  hot,  highly  conductive 
ionized  gas,  or  plasma,  streaming  outward  from  the  sun. 
This  continous  stream  of  charged  particles  la  called  the 
solar  aiafl.  The  density  of  this  'wind'  near  earth  is  on  the 
order  of  10  particles  per  cm*,  and  has  a  velocity  averaging 
300-500  km  per  second  (Jacobs).  [Ref.  6] 

Both  the  solar  wind  and  the  geomagnetic  field  exert 
pressure.  The  hot  plasma  of  the  solar  wind  pushes  against 
the  geomagnetic  field  deforming  the  field.  At  distances 
greater  than  about  13  or  14  earth  radii,  the  pressure  of  the 
solar  wind  greatly  exceeds  that  of  the  geomagnetic  field  and 
the  geomagnetic  field  will  be  swept,  along  with  the  wind. 
Prom  8  to  10  earth  radii  inward,  the  geomagnetic  field 
pressure  will  predominate,  excluding  the  solar  wind,  this 
being  the  region  of  the  J|&SL&l&21£h3£S* 

In  the  intermediate  region,  the  magnitude  of  the 
solar  wind  and  geomagnetic  field  pressures  are  comparable 
and  the  solar  wind  is  compressed  and  flows  around  the 
geomagnetic  field.  This  occurs  when  the  magnetic  energy 
density  ahead  of  the  plasma  aguals  the  kinetic  energy 
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density  of  the  streaming  plasma.  The  solar  wind  is  stopped 
at  this  point  and  forced  to  flow  around  the  magnetosphere. 
This  region  where  the  magnetosphere  starts  is  called  the 

magnetopause.  [Ref.  8] 

The  velocity  of  the  undisturbed  solar  wind  is  analo¬ 
gous  to  a  'supersonic*  velocity.  Thus  a  shock  front  is 
formed  between  the  magnetopause  and  the  solar  wind.  The 
magneto sheath  is  the  region  of  severe  turbulence  that  exists 
between  the  shock  front  and  ths  magnetopause. 

Since  the  solar  wind  always  travels  outward  from  the 
sun,  the  effect  of  the  wind  on  the  earth's  main  field  is  not 
completely  symmetric,  although  it  is  almost  symmetric  about 
an  axis  through  the  earth  and  sun.  A  geomagnetic  is 

formed  where  the  wind  sweeps  the  geomagnetic  field  along 
with  it  on  the  nightslde  of  the  earth  [Ref.  4].  Figure  1.5 
pictorially  represents  the  effects  of  the  solar  wind  on  the 
geomagnetic  field. 

5.  llAi  2t  t ‘M 

The  geomagnetic  field  changes  with  time.  As  previ¬ 
ously  mentioned,  very  slow  variations  in  the  main  field  with 
periods  on  the  order  of  years  to  thousands  of  years  Are 
referred  to  as  at SM Aft I  variati ons.  Secular  variations  are 
geologic  or  '  paleomagnetic'  in  origin.  Secular  variations 
are  not  caused  by  a  strength  or  orientation  change  of  the 
center  dipole.  Paleomagnetic  studies  are  used  to  determine 
the  secular  variation.  Geologic  structure,  especially 
conductivity  structure,  may  partially  mask  the  secular  vari¬ 
ation  at  one  point  on  the  earth  as  compared  to  that  at 
another  point. 

Other  time  variations  of  the  field  can  be  categor¬ 
ized  into  Iftliftiisa  llSl&S  and  &££££&§& 

fields.  Disturbed  variation  fields  include  gegmagnetlg 
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Figure  1.5  Configuration  of  The  Magnetosphere 


BiSES£lltiu£i£Al  which  vill  be  discussed  separately,  since 
they  are  of  particular  interest  as  a  noise  source  for  HAD 
sensors. 

a.  Quiet  Variation  Fields 

Quiet  variation  fields  ire  those  which  are  not 
due  to  disturbances  in  the  interplanetary  environment  and 
which  vary  slowly  and  regularly  [Ssf.  4]. 

There  are  several  contributing  fields  to  the 
quiet  variation.  These  include  the  Solar  Qaifit  fiflill 
(S3),  +-he  Sail!  2a£ia£i£H  (L)  *  and  the  daily 

variation  due  to  aagnetosph  eric  effects. 


The  Solar  Quiet  (Sq)  variation  is  the  name  given 
to  the  pattern  of  diurnal  field  variation  with  respect  to 
solar  local  tiae  which  is  caused  by  currents  flowing  in  the 
ionosphere  (Matsushita)  [Ref,  11].  The  aajor  portion  (about 
two-thirds)  of  the  Sq  field  is  due  to  what  is  referred  as  an 
atmospheric  drpaao.  High  speed  tidal  winds  are  generated  by 
solar  heating  causing  convection  of  the  upper  atmosphere 
(Ref.  4].  These  winds  produce  a  stationary  current  system 
by  moving  the  conducting  particles  of  the  upper  atmosphere 
across  the  geomagnetic  field  lines.  The  daily  variation  is 
caused  by  the  earth  rotating  under  the  current  system.  The 
remaining  third  of  the  sq  variation  is  caused  by  currents  in 
the  earth  induced  by  the  primary  currents  in  the  ionosphere. 

The  sq  field  can  be  shown  to  be  latitude  depen¬ 
dent  reaching  a  maximum  at  the  magnetic  equator  where  a 
concentration  of  current,  the  equatorial  electro  jet,  exists 
[Rsf.  4].  The  maximum  horizontal  component  intensity  is 
about  100  nT  at  the  equator  with  25  to  50  nT  mor9  likely  at 
higher  latitudes. 

Longitudinal,  seasonal,  and  solar  cycle  depen¬ 
dencies  also  occur  for  the  Sq  field. 

An  example  of  the  quiet-day  variation  at 
Monterey,  California  is  shown  in  Figure  1.6  and  is  summa¬ 
rized  in  Table  I,  This  data  was  taken  using  a  Cesium  Vapor 
total  field  magnetometer  in  February,  1979. 

The  Lunar  Daily  Variation,  L,  is  approximately 
one-tenth  the  magnitude  of  the  Sq  field  and  exhibits  a 
semi-diurnal  behavior  in  lunar  time  [Ref.  4].  The  major 
difference  is  that  the  winds  are  caussd  by  lunar-solar 
gravitational  tides.  The  i  field  is  dependent  on  seasonal 
influences,  lunar  phase,  the  solar  cycle,  and  latitude. 
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figure  1.6  Variation  of  Total  Geoaagnetic  Field  Intensity, 
February  26,  1  979 
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Qui«t  Day  Variation  in  La  Mesa  Village,  February,  26,  1979 
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A  diurnal  effect  due  to  the  dayside-nightside 
difference  in  coapression  by  the  solar  wind  of  the  geomag¬ 
netic  field  causes  a  small  variation  of  the  order  of  3  aT 

[  Stf  ■  4]. 

b.  Disturbed  Variation  Fields 

SiSlUCkftd  liEiiiisffl  flails  are  geomagnetic  field 
variations  that  appear  to  be  the  result  of  interplanetary 
environmental  changes  and  do  not  possess  a  simple  period¬ 
icity.  These  variations  include  ionospheric  disturbance^, 
the  asiaifl,  aafiima&is  aiaciar  ***  asaiaaaatis 
iicE2£Jiiaaii52M- 

An  ionospheric  iiatacfeaaas  is  a  departure  from 
the  normal  behavior  of  the  ionosphere. 

Sfil&E  Cl&U  a££SSta  (SEE)  are  magnetic  distur¬ 
bances  produced  by  X-rays  emitted  from  the  solar  flare. 
SFE ' s  usually  have  a  rapid  onset,  typically  a  few  minutes  in 
duration,  followed  by  a  slower  return  to  normal.  The  entire 
event  lasts  on  the  order  of  an  hour  (Reid).  [Ref.  12] 

AU£  SEAS  Ma  caused  by  the  precipitation  of 
charged  particles  down  magnetic  field  lines  into  the  atmos¬ 
phere  and  can  be  one  of  the  brightest  visual  phenomena  in 
the  3ky.  The  more  intense  and  active  auroras  occur  with 
geomagnetic  disturbances  and  greatly  increase  Ionization  as 
well  a3  creating  the  spectactular  visual  displays.  [Ref.  4] 
geomagnetic  §t2£lS  are  due  to  a  change  in  the 
dynamic  pressure  of  the  solar  wind.  A  typical  storm  begins 
with  a  compression  of  the  magnetic  field  by  an  increase  in 
solar  wind  dynamic  pressure  called  a  sadden  commencement 
<S£)  »  Sbififc  i  ncreases  the  magnetir  field  (In*  so-called 
"gradual  storm"  begins  with  a  gradual  increase  in  field 
strength).  The  increase  in  field  strength  is  on  the  order 
of  several  tens  of  nanoTeslas  (nT)  and  takes  about  one  to 
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six  minutes  to  rise.  If  a  disturbance  starts  with  an  SC  but 
lacks  the  succeeding  stages  of  a  storm  it  is  referred  to  as 
a  sudden  impulse  (§£)  .  Following  the  SC,  the  field  remains 
compressed  for  two  to  eight  hours  in  the  Phase  of 
the  storm.  The  main  follows  the  initial  phase.  Over 
a  period  of  hours  to  a  da?  a  westward  ring  currant  is  set  up 
at  a  distance  of  several  earth  radii  whose  magnetic  field 
leads  to  a  decrease  in  field  strength  on  the  order  of  100 
nT.  This  decrease  overshoots  the  agiilibrium  field  strength 
and  leads  into  the  recovery  pbasfi  of  a  da?  or  longer  where 
the  field  returns  towards  its  prestora  strength  as  the  ring 
currents  gradually  dissipate  [Ref.  4],  (Matsushita) 
fBef.  13  3.  &  magnitude-time  graph  of  a  typical  geomagnetic 
storm  is  shown  in  Pigurs  1.7. 


Figure  1.7  Typical  aidlatitmde  Geomagnetic  storm 


G «o Magnet ic  ale rommlsat ions  are  rapid  fluctua¬ 
tions  in  the  surface  magnetic  field  with  periods  of  about 
0.2  seconds  to  10  minutes  (frequencies  about  0.0016  to  5.0 
Hz)  .  These  are  observed  as  a  type  of  geomagnetic  distur¬ 
bance  by  ground  based  magnetometers  [Bef.  4]. 
Micropulsations  will  be  discussed  in  depth  later. 

6.  ElSiSJtefi  a£  39££S&£  &&&  l&ZtSL 

The  geomagnetic  field  vector  is  measjrad  or  charac¬ 
terized  at  any  point  by  its  direction  and  magnitude.  This 
can  be  done  in  terms  of  some  set  of  three  independent  param¬ 
eters  such  as  two  direction  angles  and  the  magnitude,  or 
three  perpendicular  components.  [Bef.  4] 

The  system  of  coordinates  commonly  employed  for 
describing  the  geomagnetic  field  on  the  surface  of  the  earth 
is  shown  in  Figure  1.8.  The  field  is  measured  in  terms  of. 
local  (geodetic)  coordinates  with  respect  to  True  North. 

The  various  coordinates  are  referred  to  as  magnetic 
elements  and  are  defined  as  follows: 

§:  Total  Field  Intensity  (the  symbol  £  is 
also  used) 

fl:  Horizontal  Component 

Northward,  or  North-South  Component 
J:  Eastward,  or  East-West  Component 

!J:  Downward,  or  Vertical  Component 

fi:  Declination  or  magnetic  variation. 

This  is  the  angle  between  J  and  $  and  is 
measured  positive  eastward. 

J:  Inclination  or  dip  angle.  This  is  the 
angle  between  fl  and  g  and  is  measured 
positive  downward. 


ii.  imamcim  iq  mmsis  muslin  <mj2) 

A.  DEFINITION  OF  A  MAGNETIC  ANOMALT 

A  magnetic  anomaly  is  defined  as  any  spatial  variation 
or  disturbance  in  the  geomagnetic  field  which  is  due  to 
local  causes.  Anomalies  can  be  caused  by  waves,  ore 
deposits,  sea  mounts,  and  magnetized  objects  such  as  surface 
ships  and  submarines  (Anderson)  [Bmf.  14],  For  the  purposes 
of  this  research,  magnetic  anomalies  due  to  a  simple  dipole 
field,  such  as  those  generated  by  submarines,  will  be 
regarded  as  signals,  while  other  anomalies  will  be  regarded 
as  noise  and  will  be  discussed  later. 

B.  HISTORICAL  DEVELOPMENT  OF  THE  MAGNETIC  DETECTION  OF 

SUBMARINES 

1 •  Sari?  Detection  Svstem? 

Attempts  at  finding  submerged  submarines  by  sensing 
disturbances  in  the  geomagnetic  field  date  back  as  least  as 
far  as  World  War  I.  In  19  18,  Earnest  Merritt,  at  the  Naval 
Experimental  station.  New  London  studied  the  use  of  a  fixed 
coil  type  of  detector  for  use  in  moving  boats  and  airplanes 
(Slichter).  £Bef.  15] 

MAD,  originally  known  as  Magnetic  Airborne  Detector, 
and  now  called  Magnetic  Anomaly  Detection,  had  its  begin¬ 
nings  as  an  airborne  ASW  sensor  in  late  1940  and  early  1941. 

Anomalies  in  the  geomagnetic  field  caused  by  the 
pcasence  of  submarines,  are  on  the  order  of  one  to  a  few  nT 

(gammas)  in  magnitude.  This  is  quite  small  compared  to  the 

magnitude  of  the  field  itself  (30000  -  60000  nT)  .  Two 

methods  are  generally  employed  to  measure  this  small 
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disturbance,  one  is  to  us 9  a  aia^icasiai  which  measures  the 
spatial  rate  of  change  of  the  magnetic  field  or  its 
The  second  method,  and  the  one  presently  used  by 
D.S.  Navy  aircraft  is  to  use  a  ojaqnetpm^er  to  directly 
measure  changes  in  the  magnitude  of  the  magnetic  field. 

The  British  investigated  the  use  of  a  gradiometer 
for  submarine  detection,  and  by  early  1941,  had  developed  a 
two-ccil  gradiometer  system  which  could  detect  a  submarine 
at  a  range  of  200  feet  under  favorable  conditions.  This 
range  was  considered  to  be  too  small  to  be  of  operational 
value  and  work  on  such  a  system  was  terminated  when  a  magne¬ 
tometer  system  showed  promise  (Coleman).  [Ref.  16] 

By  late  1940,  Victor  V.  Vacquier  of  the  Gulf 
Research  and  Development  Company  had  developed  a  sensitive 
satuy  qble  magnetometer  intended  for  geophysical 

(mineral)  prospecting.  The  Vacquier  magnetometer  became  the 
basis  for  further  MAD  development.  The  Airborne  Instrument 
Laboratory  of  Columbia  University  continued  the  investiga¬ 
tion  of  means  of  localizing  submerged  submarines  by  MAD. 

The  simplest  saturable  core  or  fluxgate  magnetometer 
consists  of  a  saturable  or  ferromagnetic  core  around  which  a 
coil  of  wire  is  wrapped  (see  Figure  2.1).  This  coil  carries 
a  sinusoidal  current,  I (t)  ,  which  is  large  enough  to  satu¬ 
rate  the  core  during  part  of  each  cycle.  The  inductance  of 
the  coil  will  change  as  a  function  of  the  magnetization  of 
the  cere.  The  core  magnetization,  in  turn,  depends  on  the 
instantaneous  current  in  the  coil,  and,  if  present,  the 
external  magnetic  field. 

In  the  absence  of  an  externally  applied  magnetic 
field,  magnetisation  as  a  function  of  exciting  current  is 
symmetric  around  I»0  (Figure  2.2).  An  external  magnetic 
field  parallel  to  the  core1 s  axis  will  change  the  magnetiza¬ 
tion  of  the  core  and  shift  the  magnetization  curve.  This 
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Pigure  2.1  Bl«a«nt  of  Pluxgate  Magnetometer 


change  will  cause  an  asymmetry  which  can  b9  sensed  by 
analyzing  the  harmonic  content  of  the  signal.  The  coil  can 
be  combined  with  a  stabilizer  system,  which  keeps  the 
detector  element  (coil)  aligned  with  the  geomagnetic  field. 

This  was  the  basis  of  the  MAD  magnetometers  used  in 
World  War  XI  such  as  the  AN/ASQ-1,  ASQ-1A,  and  ASQ-2 
[Ref.  16]. 

The  fluxgate  magnetometer  measures  only  the  compo¬ 
nent  of  the  external  field  parallel  to  the  axis  of  the 
ferromagnetic  core.  In  order  to  measure  the  total  field  in 
this  fashion  it  is  necessary  to  align  the  ferromagnetic  core 
along  the  earth's  magnetic  field  or  by  using  mutually 


Figure  2.2  Magnetization  Veciua  Exciting  Currant 


par pendicular  (orthogonal)  cocas.  Procision  requirements 
made  this  type  of  system  difficult  to  realize  during  World 
War  II  but  it  has  since  found  many  years  of  operational 
usage  and  has  found  widespread  use  in  geophysical  explora¬ 
tion  work  and  satellite  mapping  of  the  geomagnetic  field. 

2.  sa ill  suasa 

During  World  War  II  MAD  provided  a  passive  method  of 
detection  and  tracking  of  submerged  submarines.  Within  its 
range  (then  on  the  order  of  500  feet)  MAD  gave  a  measure  of 
surprise  to  the  attacking  aircrew.  Until  the  first  attack 
was  delivered  the  crew  of  the  submerged  submarine  might  not 
even  be  aware  of  the  aircraft  hunting  it. 

Operational  deployment  of  MAD  began  on  a  limited 
basis  in  December,  1941  with  the  installation  of  the  early 
Mark  I  MAD  in  a  blimp  at  Naval  Air  station  Lakehurst,  New 
Jersey  (Ref.  16]. 


By  the  end  of  1942,  HAD  was  operational  on  board  PBY 
Catalina  aircraft,  nicknamed  *'Madcat».n  Even  though  BAD  was 
operational  at  this  tins,  it  was  not  until  February,  1944 
that  an  initial  contact  by  MAD  led  to  the  sinking  of  a 
sub  aarine. 

During  February,  1944,  the  Hadcats  of  Patrol 
Squadron  63  were  assigned  to  fly  a  MAD  barrier  patrol  of  the 
Straits  of  Gibraltar.  On  24  February  1944,  a  Madcat  of 
vp- 6 3  detected  tl-761  by  use  of  MAD  and  commenced  tracking  to 
confirn  the  contact  as  a  moving  target.  An  attack  was 
conducted  in  conjunction  with  another  Catalina,  two 
destroyers,  and  eventually  two  other  aircraft,  and  0-761  was 
sunk.  0-392  was  sunk  after  a  similar  MAD  contact  on  16 
March,  and  on  5  May,  the  third  successful  attack  on  a  0-boat 
resulting  froa  an  initial  MAD  contact  took  place  when  0-731 
was  sunk  (OEG  Mo.  51,  Price).  [Ref.  17,  18] 

The  attacks  at  the  Straits  of  Gibraltar  demonstrated 
that  though  MAD  us  a  limited  search  rate,  there  are 
scenarios  where  it  can  be  employed  effectively  as  a  search 
sensor,  such  as  providing  a  blockade  across  a  restricted 
area  without  the  presence  of  surface  craft  (DBG  Mo.  54) 
CRef.  19]. 

3 .  £U£ia&£  2IS1SM 

The  magnetometer  system  in  current  operational  use 
is  the  optically  pumped  magnetometer.  The  optically  pumped 
magnetometer  measures  the  external  magnetic  field  by  making 
use  of  the  fact  that  when  an  atom  is  immersed  in  a  magnetic 
field,  its  energy  levels  are  split.  This  is  known  as  the 
Zeeman  effect.  For  the  fields  of  interest,  the  amount  of 
splitting  of  the  levels  is  proportional  to  the  intensity  of 
the  magnetic  field.  By  measuring  the  separation  between  the 
levels  the  magnitude  of  the  magnetic  field  can  be 
det er mined. 
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This  type  of  magnetometer  usually  makes  use  of 
Cesium  or  Rubidium  vapors,  or  Helium  gas.  The  currant  oper¬ 
ational  HAD  systea,  the  AN/ASQ-81,  is  a  Heliun  gas  optically 
pumped  magnetometer. 

In  an  optically  pumped  magnetometer  (see  Figure 
2.3),  the  sample  vapor  or  gas  (such  as  Helium)  is  collected 
into  an  absorption  cell.  circularly  polarized  light  is 
passed  through  the  cell  giving  up  some  of  its  energy  in 
exciting  or  pumping  the  electrons  of  the  sample  gas  to 
higher  energy  levels.  These  electrons  then  fluoresce  to 
lover,  metastable  energy  states.  This  is  'optical  pumping.' 
A  detector  monitors  the  degree  of  optical  pumping  by  meas¬ 
uring  the  transparency  of  the  gas  cell. 


Figure  2.3  Metastable  Helium  Magnetometer 


The  actual  separation  of  the  energy  levels  is  meas¬ 
ured  by  applying  a  weak  R.  F.  field  which  redistributes  the 


34 


electrons  aaong  the  ground  state  sublevels.  These  tran¬ 
sitions  will  only  occur  when  the  R.  ?.  field  has  a  partic¬ 
ular  frequency  (called  the  Larmor  frequency)  which  is 
proportional  to  the  separation  between  the  levels  and  there¬ 
fore  to  the  external  aagnetic  field  [Ref.  20].  Thus  aeas- 
useaent  of  this  frequency  yields  a  direct  value  of  the 
aagnetic  field. 

**  aijutti  Uftiilafiaaata 

optically  puaped  aagnetoaeters  are  the  most  sensi¬ 
tive  sensor  for  fleet  usage  today,  increasing  capability  to 
Dales  use  of  the  phenomenon  of  superconductivity  will  yield 
more  sensitive  HAD  sensors.  superconducting  Quantum 
Interference  Devices  (SQUIDS)  which  make  use  of  the 
Josephson  effect  have  a  theoretical  sensitivity  of  10-*  nT 
as  opposed  to  the  theoretical  value  of  0.001  nT  and  the 
operationally  realized  value  of  0.01  nT  in  current  fleet 
systens  (Chilton)  [Ref.  20].  Held  changes  on  the  order  of 
10-*  nT  have  actually  been  measured.  This  SQUIDS  are  used 
in  both  superconducting  magnetometers,  and  superconducting 
gra  diometers. 

A  Josephson  junction  consists  of  a  thin  layer  of 
insulator  between  two  superconductors.  This  junction  has 
the  property  that  &  current  can  flow  across  it  without 
developing  a  voltage  up  to  some  maximum  current.  A  voltage 
is  developed  for  all  currant  values  greater  than  the  maximum 
current. 

In  a  superconducting  magnetometer  a  superconducting 
ring  containing  a  pair  of  Josephson  junctions  is  used  to 
measure  the  amount  of  magnetic  flux  penetrating  the  loop. 
This  system  is  a  vector  magnetometer  which  measures  varia¬ 
tions  in  only  one  field  component. 


By  using  nor*  than  one  superconducting  loop,  it  is 
possible  to  measure  magnetic  field  gradients,  and  construct 
a  superconducting  gra disaster.  This  is  also  a  component,  or 
vector,  sensor. 

These  instruments  have  been  used  extensively  in  the 
laboratory  and  to  some  extent  in  geophysical  work,  but 
considerable  engineering  problems  remain  to  be  solved  before 
operational  Navy  use  can  be  contemplated. 

C.  NAD  SIGNAL  AND  BANDPASS 
1*  SQJ1IS*  Ql  Ills  Sisn&l 

The  NAD  signal  results  from  moving  a  magnetometer 
through  the  magnetic  field  of  a  submarine,  which  can  be 
approximated  by  the  field  of  a  magnetic  dipole.  Figure  2.4 
depicts  the  formation  of  a  submarine  caused  anomaly. 

The  magnetic  moment  of  an  object  in  the  earth's 
magnetic  field  can  be  due  to  permanent  magnetization, 
magnetization  induced  by  the  earth's  magnetic  field,  or  a 
combination  of  both.  En  the  case  of  a  submarine  hull  both 
causes  are  present  with  a  small  amount  of  permanent 
magnetization  produced  by  hull  stress  in  metal  components 
during  construction  and  stress  caused  by  submarine  diving 
and  surfacing.  This,  however,  is  a  minor  contribution.  The 
most  important  constituent  of  the  submarine  magnetic  moment 
is  due  to  magnetization  induced  in  the  hull  by  the  presence 
of  the  geomagnetic  field.  [Ref.  21] 

The  induced  field  of  a  submarine  depends  on  the 
effective  permeability  along  the  vertical,  athwartships,  and 
longitudinal  axes  of  the  submarine.  This  information  , 

taken  together  with  the  strength  and  dip  angle  of  th* 
geomagnetic  field,  yields  a  fairly  precise  calculation  of 
the  magnetic  anomaly  produced  by  the  submarine.  *  Deperming' 
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(C>  RESULTANT  DISTORTION  OF  FIELD  (d)  SIGNAL  RECEIVED  [Ref .  2(j 


of  the  submarine  hull  cancels  out  the  permanent  hull 
magnetization  leaving  the  induced  magnetization  as  t:.e 
principal  signal  source.  C  Bef.  21] 

The  HID  signal  is  approximated  as  the  projection  of 
the  submarine  dipole  field  onto  the  geomagnetic  field 
vector.  This  approximation  is  good  because  the  magnitude  of 
the  dipole  field  is  very  much  smaller  than  the  magnitude  of 
the  geomagnetic  field  (a  few  nanoleslas  as  opposed  to 
approximately  50000  &T) .  Therefore,  whenever  the  dipole 

field  is  perpendicular  to  the  earth's  field  a  region  of  zero 
signal  will  result.  The  signal  recorded  by  the  an/asq-81  or 
other  HAD  equipment  is  a  mapping  along  the  aircraft's  flight 
path  of  this  signal.  figure  2.5  qualitatively  describes 
some  aspects  of  the  MAD  signal. 

2.  Aaiama  maiiana 

The  submarine  anomaly  signal  shape  is  a  function  of 
the  dip  angle  of  the  geomagnetic  field,  the  magnetic  heading 
of  the  aircraft,  magnetic  heading  of  the  submarine  dipole, 
and  the  lateral  range  between  the  aircraft  and  the  subma¬ 
rine.  These  factors  detwermine  the  'A*  coefficients  for  the 
Anderson  functions  below. 

In  1949,  J.  E.  Anderson  of  the  Naval  Air  Development 
Center  determined  that  the  MAD  signal,  obtainad  along  any 
course,  consisted  cf  a  linear  combination  of  three  basic 
components.  Different  shaped  signals  could  by  changing  the 
proportionsl  contribution  of  these  basic  components.  The 
mathematical  representation  of  these  components  are  now 
referred  to  as  the  Anderson  functions.  Osing  the  dimension¬ 
less  parameter  • b '  (defined  as  the  distance  traveled  along 
the  aircraft  track  divided  by  the  slant  range  at  closest 
point  of  approach  (CPA)  ,  b“0  at  CPA)  the  anomaly  can  be 

represented  by 
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The  strength  of  the  detected  siqnal  is  seen  to  fall 
off  as  the  cube  of  the  distance  between  the  aircraft  and  th© 
target. 

B7  analyzing  tb.a  functions  and  their  coefficients, 
it  can  be  shown  that  the  optimum  orientation  for  maximum 
anomaly  signal  detection  is  when  the  aircraft,  target 
dipole,  and  geomagnetic  field  are  lined  up  together  as 
closely  as  possiole.  Specifically,  this  occurs  when  the 
submarine  moment  and  the  aircraft's  track  are  oriented 
North-South.  [Ref.  1ft,  21] 

The  Fourier  transform  of  the  component  functions 
were  taken  to  determine  the  frequency  distribution  of  energy 
in  the  MAD  signal.  The  anomaly  signal  components  are  shown 
in  Figures  2.6  through  2.3  (Anderson)  [Ref.  1ft],  The 
Fourier  transforms  of  the  signal  components  with  a  platform 
velocity  of  150  knots  are  shown  in  Figure  2.9. 


ill 


Figure  2.6  Basic  HAD  component  fO(b) 


Figur*  2.7  Basic  MAD  Component  fl(b) 
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3.  its  m  Eiiisi 

The  MID  signal,  being  a  transient,  is  essentially  a 
broadband  (as  opposed  to  a  discrete  frequency)  signal.  The 
frequency  spectra  of  energy  content  shows  how  the  signal  is 
contained  in  a  certain  frequency  range.  In  order  to  screen 
out  unwanted  noise,  bandpass  filtering  is  used  in  processing 
magnetic  anomaly  detections. 

The  major  factors  influencing  the  determination  of 
an  optimum  MID  filter  frequency  are  aircraft  tipeed  and  the 
slant  range  from  the  aircraft  to  the  target  at  CPA.  Factors 
exerting  a  minor  influence  include  dip  angle,  aircraft 
heading,  and  target  dipole  orientation. 

The  range. in  frequency  variation  is  about  an  octave, 
with  the  higher  values  occuring  when  the  passes  are  made 
parallel  to  the  axic  of  the  dipole  and  the  lower  when  the 
pass  is  perpendicular  to  the  dipole.  Anderson  empirically 
determined  that  the  center  frequency  of  the  filter  is  given 
by 
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where  f 
v 

and  Z 


frequency  in  Hertz 
aircraft  velocity  in  ft/sec 
greatest  anticipated  CPA 
range  in  feet 


The  ASQ-31  Bandpass  Filter  has  highpass  settings  of 
0.04,  0.06,  0.03,  and  0.1  Hz,  and  lowpass  settings  of  0.2, 
0.4,  0.6,  and  2.0  Hz  (Orion  Service  Digest  26)  [Bef.  22]. 
There  are  no  recommended  settings  for  normal  operation  since 
background  noise  varies.  The  filter  characteristics  for  the 
0.06  to  0.6  Hz  settings  are  shown  in  Figure  2.10  (Orion 
Service  Digest  28)  [Ref.  23], 
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Tha  figure  is  a  representation  of  the  adjustable  band  pass 
filter  of  the  ASQ-81  by  itself.  Other  parts  of  the  ASQ-3'1 
system  add  in  an  additional  high  pass  filter  which  adds 
another  12  dB/octave  roll  off  to  the  low  frequency  end  in 
figure  2.10.  Thus  Figure  2.10  would  represent  the  charac- 
teristics  of  the  entire  ASQ-81  syste®  if  the  roll  off  of  36 
dB/octave  is  changed  to  46  dB/octave.  [Ref.  24] 


Figure  2.10  ASQ-81  Filter  Characteristics:  0.06  -  0.6  Hz 

Using  equation  2-2,  the  center  frequency  for  a  MAD 
signal  with  a  CPA  slant  range  of  300  feet  and  aircraft 
velocity  of  220  knots  is  0.49  Hz.  Using  a  lowpass  filter 
setting  above  0.6  Hz  would  not  be  very  useful. 

One  reason  tha  ASQ-81  filter  extends  up  to  2.0  Hz 
instead  of  just  0.6  Hz  is  the  following.  The  early  Service 
Test  Engineering  Model  (STEM)  ASQ-81  that  preceeded  the 


currant  ASQ-81  production  model  had  two  bandpass  filters, 
one  adjustable  from  .04  to  .6  Hs,  and  the  other  fixed  at  .75 
to  10  Hz.  During  testa,  the  .75  to  10  Hz  channel  proved  to 
be  helpful  in  monitoring  the  STEM  A8Q-31  "system"  noise 
level.  During  design  reviews  for  tho  production  ASQ-31,  it 
was  decided  that  it  night  be  useful  to  retain  the  ability  to 
monitor  "system"  noise.  In  order  to  do  this  easily,  the  .75 
to  10  Hz  band  was  dropped  and  a  2  Hz  position  was  added  to 
the  adjustable  filter.  [8ef.  24] 

The  center  frequency  at  1200  feet  and  160  knots  is 
.09  Hz.  Consequently,  the  highpass  filter  settings  could 
affect  slower,  longer  range  HAD  detections. 

The  choice  of  the  highpass  (lover  end  of  bandpass) 
settings  have  to  be  made  after  considering  the  noise  present 
at  time  of  system  operation. 
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saaasss  at  me  naias 


A.  INTRODUCTION 

On  a  practical  level  MAD  noise  la  defined  aa  magnetic 
disturbances  falling  within  the  MAD  pasaband  (0.04  to  2.0 
Hz,  or  sines  0.6  Hz  is  tbs  normal  uppsr  limit,  0.04  to  0.6 
Hz)  and  hawing  an  amplitude  grsatsr  than  0.01  n&noTeslas 
(the  system  sensitivity  of  the  AN/ASQ-31  MAD  system)  . 

MAD  noise  sources  can  be  divided  into  the  following 
categories : 

-Equipment  Noise 
-Aircraft  Platform  Noise 
-Aircraft  Manuever  Noise 

-Gradient  Noise  Due  to  Aircraft  Motion  Through  the 
Geomagnetic  Field 
-Geologic  Noise 

-Noise  from  find  Naves  and  Swells 
-Geomagnetic  Noise 

B.  SENSOR,  PLATFORM  AND  MANUEVER  NOISE 

Sensor  noise  is  the  self-noise  generated  by  the  opera¬ 
tion  of  the  equipment  itself.  This  can  be  partially  due  to 
the  fact  that  the  detector  element  is  misaligned  with 
respect  to  the  geomagnetic  field  vector.  Changes  in  lamp 
intensity,  photodetector  noise  and  noise  in  the  electronic 
circuits  can  also  contribute  to  sensor  noise  [Ref.  21].  The 
self-noise  limitation  of  the  ASQ-81  is  0.01  nT. 

Platform  noisa  is  generated  by  components  fixed  to  the 
aircraft  in  the  vicinity  of  the  sensor  [Ref.  23]. 
Permanent,  induced  and  eddy-current  magnetic  fields  ars» 
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associated  with  the  airframe.  Permanent  magnetic  fields  are 
due  to  aircraft  structure  or  equipment  having  ferromagnetic 
parts.  This  field  changes  its  orientation  with  respect  to 
the  geomagnetic  field  vector  as  the  aircraft  manuevers, 
causing  field  fluctuations  near  the  magnetometer. 

Platform  noise  is  also  induced  in  aircraft  ferromagnetic 
structures  by  the  geomagnetic  field.  Similarly,  eddy-cur¬ 
rants  are  induced  in  aircraft  skin,  ribs,  and  frames,  and 
these  currents  ,  in  turn,  cause  additional  magnetic  fields. 
Thus,  rapid  aircraft  v manuevers  will  Induce  changes  in  the 
magnetic  field  sensed  by  the  magnetometer. 

Platform  noise  in  aircraft  mounted  sensors  is  countered 
by  applying  equal  and  opposite  magnetic  fields  to  the  sensor 
in  a  process  called  compensation.  Towed  MAD  systems  are 
essentially  frre  from  this  type  of  noise. 

C.  GRADIENT  NOISE 

Gradient  noise  can  be  divided  into  turn  noise  and  noise 
due  to  changes  in  altitude. 

Turn  noise  is  a  problem  when  ' HAD  trapping*  or  using  MAD 
for  tracking  a  target.  The  earth's  magnetic  field  has  a 
horizontal  gradient  (in  this  case  the  magnitude  varying  with 
latitude).  As  the  aircraft  moves  in  the  direction  of  the 
gradient  the  field  strength  changes.  The  r.oise  due  to  a  MAD 
trapping  or  hunting  circle  is  centered  in  frequency  at  the 
reciprocal  of  the  time  taken  to  complete  one  revolution.  In 
the  case  of  a  two-ainute  circle,  the  noise  would  be  centered 
at  0.00333  Hz,  well  below  the  filter  used  in  the  fiSQ-81. 
The  horizontal  gradient  noise  due  to  flying  a  cloverleaf 
pattern,  for  the  most  part,  also  falls  below  the  MAD 
paasband.  [Ref,  21] 

Vertical  gradient  noise  is  due  to  changes  in  sensor 
altitude.  An  altitude  gradient  of  up  to  0.005  to  0.01  nT 


per  foot  exists  in  the  earth's  main  field.  la  areas  of 
geological  anomalies  this  gradient  is  even  larger.  Fast 
altitude  or  aircraft  pitch  changes  can  cause  a  magnetic 
field  fluctuation  of  sufficient  amplitude  to  be  of  concern. 
To  avoid  vertical  gradient  noise,  altitude  compensation 
equipment  is  used.  [Ref.  23] 

D.  ENVIRONMENTAL  NOISE 

Magnetic  noise  from  sources  existing  in  the  natural 
environment  include  geologic  noise,  temporal  variation  in 
the  earth's  magnetic  field,  and  noise  due  to  ocean  waves  and 
swells. 

Geologic  noise  has  its  source  in  naturally  occuring 
magnetic  anomalies  caused  by  magnetic  material  present  in 
the  earth's  crust.  when  the  sensor  passes  over  geological 
anomalies,  the  relative  motion  causes  a  MAD-liks  signal  to 
be  recorded.  Geologic  noise  is  usually  more  pronounced  in 
shallow  water  as  the  sensor  is  much  closer  to  the  source  of 
the  noise.  Geological  magnetic  anomalies  are  often  associ¬ 
ated  with  such  oceanographic  features  as  seamounts  and  ocean 
ridges. 

Se a  water  is  a  conducting  medium  which  is  transported  by 
the  physical  motion  of  water  wavra  la  the  presence  of  the 
geomagnetic  field.  This  motion  induoos  currents  in  the  sea. 
These  currents  give  rise  to  secondary  magnetic  fields,  which 
add  vectorially  to  the  quasiatatic,  geomagnetic  field 
(Weaver)  [R*f„  25].  These  fields  can  be  detected  at  signif¬ 
icant  distances  above  the  sea  surface  and  fall  off  exponen¬ 
tially  with  altitude.  Figure  3. 1  is  a  plot  for  several 
surface  wave  periods  of  the  Induced  magnetic  field  per  meter 
amplitude  of  tho  surface  wave.,  These  induced  fieldn  can  be 
a  problem  at  the  low  altitudes  where  MAD  is  U3*d. 


Figure  3.1  induced  Magnetic  Field  per  Beter  of  iavaUeight 
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S.  GEOMAGNETIC  NOISE 


Temporal  variations  in  the  earth's  magnetic  field  vith 
frequencies  in  the  MAD  bandpass  and  amplitudes  greater  than 
0.01  nT  have  become  known  as  SlfiliS&aiifi  noise  in  the  MID 
literature. 

Quiet  daily  variations,  such  as  the  sq  and  L  variations, 
have  periods  sufficiently  long  to  fall  far  below  the  MAD 
passband. 

Geomagnetic  storms  have  been  discussed  previously. 
Bap  id  fluctuations  vith  high  amplitude  and  falling  within 
the  MAD  passband  occur  in  connection  with  geomagnetic 
storms. 

Geomaqnat ic  micro pulsations  comprise  the  last  category 
of  geomagnetic  noise  to  be  discussed. 

1*  SiSMOUJtUS 

Geomagnetic  micropulsat ions  are  rapid  fluctuations 
of  the  earth's  magnetic  field  with  periods  from  0.2  seconds 
to  10  minutes  and  amplitudes  from  about  0.1  nT  to  as  high  as 
a  few  tens  of  nT's.  These  fluctuations  are  caused  by  elec¬ 
tromagnetic  perturbations  propagating  in  the  magnetosphere 
as  hydro  magnetic  waves  (Nlshida).  [Ref.  26) 

Micropulsations  are  classified  by  morphology,  that 
is,  by  examining  periods,  amplitudes,  times  of  occurence  and 
other  observed  characteristics.  Micropulsations  have  been 
placed  into  two  broad  categories:  i£E£SLJllAl  and  continuous. 
The  irregular  pulsations  are  represented  by  the  symbol  Pi, 
and  the  continuous  pulsations  by  Pc.  Table  II  shows  the 
breakdown  by  period  of  the  various  micropulsation  classes. 

It  can  be  seen  that  for  MAD  geomagnetic  noise  the 
lower  frequency  Pci,  Pc2,  Pc3,  and  Pil  pulsations  are  of 
interest. 
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TABLE  II 


Geomagnetic  Micropulsation  Classes 


NOTATION 
Pci 
PC  2 
Pc  3 
Pc  4 
Pc  5 
Pil 
Pi  2 


PERIOD  (sac) 


AVERAGE  AMPLITUDE  (n?) 
0,05  -  0.1 
0.1  -  1.0 
0.1  -  1.0 
0.1  -  1.0 
1  -  10 
0.05  -  jj.1 

(after  Jacobs) 


pci  pulsations  are  regular  sinusoidal  oscillations 
with  periods  normallly  falling  in  the  0.3  to  4  second  (0.25 
to  3.33  Hz  frequency)  range.  Pol's  may  start  as  separate 
bursts  and  gradually  develop  into  a  series  of  pulsations 
which  could  last  for  hours.  They  may  also  occur  as  consecu¬ 
tive  groups  of  pulsations  with  sharply  varying  frequencies. 
The  average  amplitude  of  Pci's  is  0.05  to  0.1  nT  and  they 
tend  to  have  a  single  well  defined  frequency.  Pci  pulsa¬ 
tions  with  frequency  less  than  0.5  Hz  are  more  common  at 
high  latitudes  than  at  aid  latitudes.  These  pulsations 
occur  in  the  daytime  in  the  auroral  zone,  and  at  night  and 
early  morning  hours  in  lower  latitudes.  Pci's  are  charac¬ 
teristic  of  the  quiet  and  weakly  disturbed  states  of  the 
geomagnetic  field  and  show  an  increase  in  activity  one  to 
two  hours  before,  and  four  to  seven  days  after,  a  magnetic 
storm.  (Ref.  6] 

Pc2  and  Pc3  pulsations  are  grouped  together  in  most 
characteristics.  The  amplitudes  of  the  oscillations  are 
usually  under  0.5  nT  and  the  typical  frequency  range  is  0.03 
to  0.2  Hz.  These  pulsations  are  normally  found  during  the 
day  'with  their  activity  reaching  a  maximum  around  noon.  Pc3 
pulsations  show  a  seasonal  variation  with  a  minimum  of 
activity  occuring  during  winter. 


Pc4  and  Pc5  are  large  amplitude  fluctuations,  bat 
fall  below  the  frequency  band  of  interest  for  MAD. 

Pil  pulsations  have  ar.  irregular  form  with  an 
average  amplitude  of  0.01  to  0.1  nT  and  a  frequency  mainly 
in  the  0.10  to  0,17  bz  range.  spectral  analysis  of  these 
pulsations  show  a  braod  band  of  frequencies.  Pil  aaplitudes 
have  maximum  values  in  the  auroral  zones,  with  the  intensity 
of  the  pulsations  decreasing  vith  decreasing  latitude. 
Pil's  are  normally  observed  in  the  late  night  and  early 
morning  hours,  and  show  an  increase  in  activity  with 
increased  geomagnetic  field  disturbance.  £H»f„  6] 

Geomagnetic  micropulsations  can  be  observed  anywhere 
on  the  globe,  at  various  times  of  day  and  year,  and  in  both 
quiet  and  disturbed  geomagnetic  field  conditions. 
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iv.  2E  H&UIA1X8S  £S25i  ASCIIS  hZlIUU 

A.  INTRODUCTION 

Previous  chapters  have  defined  geomagnetic  activity  as 
it  applies  to  Magnetic  Anomaly  Detection.  In  this  chapter, 
methods  of  evaluating  that  activity  will  be  examined, 
including  methods  currently  in  use  by  the  fleet. 


B.  GEOMAGNETIC  INDICES 

A  geomagnetic  index  is  sinply  a  measure  used  to  quantify 
and  describe  time  variations  of  the  earth's  magnetic  field 
resulting  from  solar-terrestrial  relationships.  These 
indices  are  commonly  used  to  express?  the  intensity  and 
depict  the  character  of  geomagnetic  activity  throughout  the 
day. 

For  the  most  part,  geomagnetic  indices  developed  as 
range  indices,  measuring  the  difference  between  the  high  and 
low  values  for  different  field  components  measured  during 
the  day  by  magnetic  observatories  (Lincoln)  [Ref.  27].  Most 
current  indices  are  of  the  range  type,  but  other  indices 
have  been  developed  which  are  more  subjective  or  qualitative 
in  nature. 

Geomacrnetic  indices  are  designated  by  a  letter  coda  such 
as:  C,  Ci,  Cp,  C9 ,  Q,  ft,  W  measure,  Dsx,  K ,  its,  Kp,  ak,  Ak, 
ap,  and  Ap.  There  are  add! tonal  indices  in  use. 

The  C,  Ci,  Cp,  and  C9  indices  are  daily  magnetic  field 
character  figures.  The  C  index  is  the  daily  character 
figure  for  a  single  observatory.  In  this  scale,  C*0  indi¬ 
cates  a  quiet  day,  c* 1  a  moderately  disturbed  day,  and  C*2  a 
heavily  disturbed  day.  The  daily  international  character 
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figure,  Ci,  is  the  arithmetic  mean  of  the  c  indices  reported 
by  participating  observatories  around  the  world.  Cp,  the  ;| 

daily  planetary  character  figure,  is  similar  to  Ci,  except  .  , 

that  it  is  derived  free  the  values  of  Kp  and  ap.  C9  is  a 
contracted  scale  for  Ci  an  cp  with  single  digit  values 

■  ■  .•< 

running  frost  0  to  9  (Bartels).  [Hsf.  28 ]• 

The  o  and  B  indices  are  quarter- hourly  and  hourly  range  j 

indicts  respectively,  taken  at  high  latitude  stations  only. 

The  w  measure  is  an  index  of  the  equatorial  electro jet.  ; 

Dst  is  a  measure  of  ring  current  effect.  They  are  both  | 

aaplitude  indices,  (fief.  27] 

I*  £*  au  aai  &  loiisss 

The  K,  Ks,  Kp,  ah,  Ak,  ap,  and  Ap  indices  comprise  a 
group  of  related  3-hour  range  indices.  The  K  index  is  a 
single  station  coda  using  a  quasi-logarithaic  scale  front  9 
to  9  to  measure  geomagnetic  activity.  The  value  of  K  is 
determined  by  first  determining  the  difference  between  the 
lowest  and  highest  deviations  from  the  regular  daily  varia¬ 
tion  (sg)  during  a  3-hour  period.  This  range  (in  nT)  is  j 

converted  to  the  K  scale  based  on  bhe  historical  activity 
ranges  at  the  particular  observatory  involved  [Bef.  27]. 

The  conversion  for  the  Fredericksburg,  Virginia  observatory 

is  given  in  Table  III.  This  conversion  can  also  be  applied  ] 

no  the  USAF/NOAA  observatory  in  Boulder,  Colorado  C8«i*  29].  ! 

j 

The  Ks  index  is  a  standardized  K  index  which  is  i 

1 

freed  from  local  variations  and  is  then  used  to  determine 
the  planetary  3-hour  index,  Kp. 

The  equivalent  three-hour-range,  ak,  is  a  conversion 
of  the  K  index  as  shown  in  Table  IV.  In  order  to  detemine 
the  units  of  ak  for  a  particular  observatory,  divide  the 
lower  range  limit  of  K»9  by  250.  Thus  for  Fredericksburg 
and  Boulder,  ak  is  in  2-nT  units. 
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TABLE  III 


Conversion  fro  a  Bangs  to  R  for  Pr®derickaburg,  Va, 


K 

0 

1 

I 


8 

9 


nge^nT) 


(after  Lincoln) 


TABLE  IV 

Equivalent  Range  ak  foe  Given  K 


K  0  1  2  3  456789 

ak  0  3  7  15  27  48  80  .140  240  400 

(after  Lincoln) 


Ak  is  the  equivalent  daily  amplitude  and  is  the 
average  of  the  eight  daily  ak  values  at  a  particular  observ¬ 
atory.  This  index  is  promulgated  using  the  name  of  the 
observatory,  the  Ak  index  for  Fredericksburg  is  known  as  the 
A-Fredericksburg  or  A-Fred  index. 

The  equivalent  planetary  amplitude,  ap,  is  deter¬ 
mined  from  the  Kp  index  in  a  fashion  similar  to  that  of 
determining  ak  from  the  K  indices.  The  eight  ap  values  for 
a  given  day  can  then  be  averaged  into  the  daily  equivalent 
planetary  amplitude  Ap.  These  two  indices  are  given  in  2-nT 
units. 
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C.  GEOMAGNETIC  INDICES  IN  FLEET  USB 


1*  £il££fl]&  US&SS 

Geomagnetic  indicts  of  interest  in  connection  with 
NAD  operations  are  the  K,  Ak,  and  Ap  indices. 

fleet  operators  utilize  the  "alpha  index*  for 
per  dieting  geosagnetic  activity  over  the  entire  world 
( Ref*  30].  This  index  is  promulagated  by  the  Fleet 
Numerical  Octnaographic  Canter,  Monterey,  California,  in  the 
environmental  briefings  received  by  aircrew  personnel.  This 
index  is  the  Ap  index  as  sent  out  from  the  Space 

Environmental  Services  Canter,  Boulder,  Colorado  in  the 
Joint  USAP/NOAA  Report  of  Solar  and  Geophysical  Activity 
{Ref.  31]. 

The  Boulder  K  index  is  available  to  interested 

parties  by  telephone  recording  and  in  the  WWV  and  WWVH  radio 

broadcasts  and  is  therefore  available  to  fleet  users 

(Ref.  32]. 

The  Ak  index  for  the  Fredericksburg,  Virginia 
observatory  hna  been  used  in  studies  of  geomagnetic  activity 
as  applied  to  Magnetic  Anomaly  Detection  (Ref.  33  ]„ 

2  •  2!kft2E&iSAl  2l  &  ASd  &  I&&e££I  &  M2 

With  the  widest  useful  filter  settings,  the  MAD 
bandpass  ranges  from  0.04  to  0.6  Hz  (1.7  to  25  seconds  in 
period).  Aa  such,  in  order  for  a  geomagnetic  index  to  be 
directly  applicable  for  MAD  use,  it  should  be  sensitive  to 
that  frequency  range. 

The  K  indices,  and  the  K-deriued  i  indices,  are  not 
especially  sensitive  to  the  MAD  range.  Mayaud  [Ref.  34] 
indicates  that  these  indices  are  mainly  sensitive  to  fluctu¬ 
ations  whose  periods  are  much  longer  than  the  lower  end  of 
the  frequency  ge  analyzed,  that  is,  a  frequency  corre¬ 
sponding  to  a  pe*  .»d  of  45  minutes  (0.0004  Hz). 


On «  reason  for  this  lack  of  sensitivity  for  MID 
bandpass  geomagnetic  noise  is  that  the  amplitude  of  geomag¬ 
netic  fluctuations  varies  inversely  with  frequency,  so  that 
the  amplitude  of  the  fluctuation  increases  as  the  as  the 
frequency  decreases.  It  can  therefore  be  seen  that  the 
fluctuations  with  periods  of  an  hour  or  greater  largely 
determine  the  variation  range  used  to  calculate  the  K  index. 
The  activity  driving  the  K  and  A  indices,  will,  because  of 
band  pass  filtering,  not  even  be  observed  by  the  MAD  system, 
and  the  activity  of  interest  to  MAD  might  not  influence  the 
K  or  A  indices  at  all.  It  can  be  concluded  that  there  is  no 
direct  physical  link  between  MAD  geomagnetic  noi3e  and 
either  the  K  or  A  indices.  (Ref.  35] 

3.  2l  1  t  £0&£££  iiiil  &A2 
£UL&  HaiSft 

a.  ASQ-10A  Study 

Brennan  and  Salts  CRef.  33]  found  that  geomag¬ 
netic  micropulsation  activity  was  recorded  at  their  ASQ-10A 
MAD  magnetometer  site  in  Maryland,  when  the  A-Predericksburg 
index  was  greater  than  25.  This  occured  everytime  they  were 
recording  data  with  A-Fred  greater  than  25.  It  is  important 
to  note  that  they  observed  additional  activity  during  some 
periods  when  A-Fred  was  less  than  25. 

The  Brennan  and  Smits  study  tends  to  validate 
the  use  of  the  A  indices  as  at  least  qualitative  indications 
of  geomagnetic  noise  in  the  MAD  bandpass.  Their  study, 
however,  was  specific  to  the  ASQ-10A  magnetometer  system 
which  has  a  sensitivity  of  0.1  nT  as  opposed  to  the  0.01  nT 
sensitivity  of  the  ASQ-81  system.  The  effect  of  the  ASQ-10 
sensitivity  is  to  filter  out  most  Pci  pulsations.  ks  Pci 
pulsations  do  not  correlate  well  with  the  A  and  K  indices, 
the  filtering  out  of  these  pulsations  would  tend  to  increase 
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the  reliability  of  the  A-Fred  or  &p  indices  as  measures  of 
MAD  geomagnetic  noise.  Mason  [Bef.  36]  stated  that  the 
"occurence  of  Pci  is  well  known  to  be  associated  with  low  Kp 
values."  It  should  be  noted  that  an  operational  drawback  of 
the  filtering  out  of  Pci  pulsations  by  the  ASQ-10A  is  that 
the  system  also  filters  out  valid  Signals  of  less  than  0.1 
nT  amplitude. 

While  information  has  been  presented  that 
suggests  the  the  A-Fred  index  can  be  useful  for  MAD  geomag¬ 
netic  noise  evaluation  for  the  AN/AS3-1QA  system,  sufficient 
data  was  not  presented  to  draw  conclusions  for  index  usage 
with  the  AN/ASQ-81  system. 

b.  1SQ-81  Study 

For  two  weeks  in  April  1976,  Naval  Air 
Development  Center  personnel  operated  a  geomagnetic  observa¬ 
tory  at  the  Atlantic  Undersea  Test  and  Evaluation  Center  in 
the  Bahamas.  The  primary  magnetometer  used  for  this  observ¬ 
atory  was  the  ASQ-31  magnetometer.  One  purpose  of  this 
observatory  was  to  compare  the  K  index  (as  determined  by  the 
San  Juan  Observatory)  with  geomagnetic  activity  in  the  MAD 
band  pass.  The  conclusion  of  this  study  was  that  the  K-San 
Juan  index  did  not  correlate  with  geomagnetic  noise  in  the 
MAD  band  pass.  [Bef.  37] 

We  have  undertaken  a  correlation  analysis  using 
the  Fisher  Z  transformation  [Bef.  38]  of  the  NADC  data 
furnished  to  us  by  Ochadlick.  The  data  used  is  listed  in 
Tables  V  and  VI. 
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TABLE  7 

Obaarvad  1UTEC  Data  and  Indicaa,  April  1 1- 1 8 #  1976. 


Data 


Tia«IZ) 


m 

II 


m 

a 

03  00 
06  00 
0900 
1200 
1500 
18  QQ 


09  00 
1200 
1500 
2100 


K-S&n 

1 


Juan  K 


Maximus 

iplituda 


06 

3:84 

0.13 


0.26 

0.16 

Q.Q7 


7  y«  i  j 

4  0.19 

4  0.30 

8  0.11 

8  0.48 

I  8:8? 

8  0.19 

1  ;il 

5  0.11 

5  0.  14 

5  0.10 

5  0.04 

5  0.08 

5  0.  10 

!  8:88 

2  0.18 

2  0.10 

2  0.10 

2  0.07 

(after  Ochadliclc) 


1  ."-v.  «■ 


TABLE  VI 


observed  AO TEC  Data  and  Indices,  April  19-24,  1976, 


Pate 

19 


*i  MiS) 


K-San  Juac 


Baximum 
Amplitude  (nT) 


50  U.  It 

(after  Ochadlick) 


The  largest  ptak-to-peak  fluctuations  observed 
on  the  ASQ-81  magnetometer  in  a  three  hour  period  was 
compared  to  the  K-San  Juan  index  for  that  period.  Data 
points  with  fluctuations  greater  than  1.0  nT  were  deleted  as 
the  observatory  did  not  record  that  information.  The  sample 
coefficient  of  linear  correlation  was  0.34  with  the  .95 
confidence  interval  for  the  actual  correlation  coefficient 


being  from  0.14  to  to  0.52.  Sample  size  wee  86.  This  indi¬ 
cates  that  there  is  at  moat  a  weak  correlation  between  the 
observed  data  and  the  K-San  Juan  index.  i  much  greater 
correlation  would  bo  required  for  the  K  index  to  be  of  any 
significant  value  for  MAD  operational  use. 

The  coefficient  of  correlation  has  values  from 
-1.0  to  1.0.  A  value  of  -1.0  or  1.0  indicates  perfect 
negative  or  positive  correlation,  respectively.  A  value  of 
zero  signifies  no  .correlation  at  all. 

A  siailar  correlation  analysis  against  tho  Ap 
index  was  conducted.  The  sample  correlation  coefficient  was 
0.51  with  the  .95  confidence  interval  for  the  actual  coeffi¬ 
cient  of  correlation  being  from  0.32  to  0.65.  Again,  this 
signifies  that  only  a  weak  correlation  was  observed. 

The  weak  correlation  between  observed  geomag¬ 
netic  activity  and  the  K-san  Juan  index  suggests  that  this 
index  would  not  be*  very  useful  in  descriLing  MAD  geomagnetic 
noise,  as  this  K  index  reflected  activity  similar  to  that  in 
the  band  of  interest  only  about  one-third  of  the  time.  The 
stonger  correlation  of  the  Ap  index  indicates  that  it 
reflected  activity  similar  to  MAD  geomagnetic  noise  about 
one-half  of  the  time.  This  is  still  not  a  very  good  indica¬ 
tion  of  what  is  going  on  in  the  MAD  band  pass. 

c.  Power  Spectral  Density  Evaluation 

As  part  of  ongoing  research  at  the  Naval 
Postgraduate  School,  geomagnetic  activity  data  in  the  range 
of  the  MAD  band  pass  has  been  collected  and  analyzed  in  the 
form  of  power  spectral  density  (PSD)  curves, 

A  measure  of  MAD  band  activity  has  been  devel¬ 
oped  by  integrating  under  the  PSD  curve  and  using  compromise 
band  pass  limits  of  0.05  and  1.0  Hz.  This  type  of  index  is 
discussed  in  depth  later. 
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Two  sets  of  data  war*  analysed.  on*  sat  was 
taken  using  an  induction  coil  to  neasura  the  fluctuations  in 
one  direction  during  day,  I960  [Rtf.  39]  and  August  to 
October,  1981  [Ref.  40].  The  other  data  set  was  collected 
using  a  Cesium  vapor  toatai  field  eagnetoseter  during  the 
July  to  October,  1980  period  CRef.  41  ]. 

Correlation  analysis  of  the  single-coil  systea 
data  yielded  a  sample  correlation  coefficient  of  0.40  1  for 
the  K-Fredericksburg  index,  0.180  for  the  A-Prederieksburg 
index  and  0.155  for  the  Ap  index.  Sanple  size  ttaa  9*  .'she 
.95  confidence  intervals  for  the  correlation  coefficient  for 
K-?red,  A-Fred,  and  Ap  were  -.36  to  .85,  -.55  to  .76,  and 
-.56  to  .75,  respectively.  The  lata  for  this  test  iu 
presented  in  Table  ?! I. 


TABLE  VII 


Single  Coil  RMS  Hoist  Data  and  Indices 


Date 


A 


]} 

17 

7 

7 

7 


The  Cesium  vapor  magnetometer  data  yielded 
sample  coefficients  of  correlation  of  0.552,  0.374,  and 

0.444  for  the  K-Fred,  A-?red,  and  Ap  indices,  respectively. 
The  sample  sire  was  14.  The  .95  confidence  Intervals  were 
for  K-Fred  from  .03  to  .84,  A-Pred  from  -.20  to  .76,  and  for 
Ap  from  -.12  to  .79.  The  Cs  vapor  magnetometer  data  is 
presented  in  Table  ?l II. 
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TABLE  VIII 

C«  Vapor  RMS  Hois*  Data  and  indices,  Jul-0ct,1980 


Although  the  sample  alia 9  used  wart  too  small  to 
draw  any  meaningful  conclusions,  therm  is  lit, tit  evidence  to 
suggest  that  any  of  the  K-Pred,  A-Fr«d,  or  Ap  indices  Is  a 
very  accurate  measure  of  geomagnetic  noise  in  the  ASQ-81  mad 
band  pass. 

i 

d.  Correlation  conclusions 

Although  some  weak  correlation  does  exist 
between  the  K  indicia,  the  A-Fred  index,  the  Ap  index  and 
geomagnetic  noise  in  the  MAD  band  pass,  this  correlation  is 
incidental  and  indirect,  being  thu  result  of  a  correlation 
between  the  activity  in  the  MAD  band  and  in  the  lower 
frequency  activity  that  influences  the  K  and  A  indices. 
These  indices  are  not  directly  influenced  by  activity  in  the 
MAD  band.  The  correlation  that  does  exist  does  not  appear  to 
be  sufficiently  high  to  enable  these  indices  to  yield  accu¬ 
rate  indications  of  the  actual  MAD  band  activity.  The  use 
of  these  indices  for  anything  except  the  roughest  qualita¬ 
tive  estimation  of  active  in  the  MAD  band  pass  is  not 
recommended. 
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0.  P  BOP  OS  2D  GEOMAGNETIC  INDICES  FOR  HAD 

Overall  geomagnetic  acitivity  is  analyzed  in  both  the 
tine  and  frequency  domains.  Geomagnetic  noise  indices  for 
the  HAD  band  pass  could  be  developed  in  either  of  these 
domains. 

The  spatial  coherence  of  HAD  geomagnetic  noise  has  not 
yet  been  adequately  determined.  This  information  would  be 
necessary  in  order  to  determine  the  number  and  location  of 
mini-observatories  for  an  operational  had  noise  index. 

1*  liflft  SftElAft  AAfliXfili 

one  way  to  develop  an  index  of  had  geomagnetic 
activity  would  be  to  establish  mini-observatories  near  bases 
from  which  MAD  operations  are  conducted.  These  observato¬ 
ries  would  use  ASQ-81  magnetometers  or  different  magnetome¬ 
ters  with  ASQ-81  filter  networks,  and  could  in  real  time 
record  the  geomagnetic  noise  in  the  HAD  band.  A  measure 
such  as  the  maximum  peak-to-peak  (or  possibly  the  average 
peak-to-peak)  noise  in  a  given  time  period  could  then  be 
disseminated  to  flight  crews  operating  in  the  area  covered 
by  that  index,  obviously,  the  spatial  coherence  of  MAD  band 
activity  is  important  in  making  such  a  system  work.  This 
type  of  mini-observatory  has  been  suggested  by  References  31 
and  35. 

2.  iiaamsz  Osiaia 

Present  fleet  procedures  examine  HAD  noise  such  as 
system  and  manuever  noise  in  terms  of  the  amplitude  of  the 
fluctuation  [Ref.  30].  An  index  of  geomagnetic  noise  in  the 
HAD  band  pass  would  therefore  be  of  greatest  usefulness  to 
the  fleet  operator  if  it  were  in  units  of  the  amplitude  of 
the  signal  as  seen  by  the  HAD  equipment. 
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The  method  proposed  to  derive  a  MAD  noise  index  in 
the  frequency  domain  begins  with  obtaining  the  power  spec¬ 
tral  density  of  the  activity  in  the  SAD  band  by  Fourier 
analysis  of  the  tine  series  data  input  from  the 
magnetometer. 

By  inter grating  under  the  PSD  curve  over  the  limits 
of  the  MAD  bandpass,  a  value  in  units  of  amplitude2  will 
result.  Taking  the  square  root  of  this  value  will  yield  an 
HMS  amplitude.  Equation  4-1  represents  the  derivation  of 
this  index. 


MAD  index 


(r 


PSD  (f)  df 


) 


1/2 


where  NAD 
and  PSD(f)1 


Index  is  in  nT  (gammas) 

-  frequency  (Hz) 

■  upper  bandpass  limit 

*  lower  bandpass  limit 

*  power,  spectral  density 
(nTVhzf 


(4-1) 


The  characteristics  of  the  HAD  filter  (the  filter 
not  being  an  ideal  bandpass  filter)  could  be  applied  prior 
to  the  integration.  The  integration  it  self* could  be  done  by 
either  a  point -by- point  numerical  integration  or  by  first 
modelling  the  PSD  curve  by  polynomial  curve  fitting  and  then 
integrating  the  polynomial  over  the  range  of  the  band  pass. 
It  is  anticipated  that  this  could  be  done  in  close  to  real 
time  by  a  digital  computer,  possibly  by  a  desk  top  computer 
such  as  the  HP9845. 

The  type  of  sensor  utilized  could  be  the  ASQ-81 
magnetometer,  other  total  field  magnetometer,  or  possibly  an 
orthogonal  3-coil  system  whose  signals  can  be  combined  to 
yield  the  projection  on  the  total  field  vector  of  the  fluc¬ 
tuations.  A  single  coil  system  oriented  in  the  direction  of 
the  earth's  magnetic  field  vector  could  also  be  used. 
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A  three  coil  s ysta®  which  is  used  to  yield  an  RMS 
amplitude  is  currently  in  research  use  at  the  Naval 
Postgraduate  School, 

3.  £L£3i£tiSA5  £|  S£2Siaaa£iA£  to&ixiil 

The  proposed  indices  discussed  above  are  intended  to 
be  real  time  aeasures  of  the  geomagnetic  noise  in  the  MAD 
bandpass.  Whether  or  not  such  activity  can  be  predicted 
ahead  of  time  needs  to  be  looked  into. 

While  there  is  no  aodel  for  the  background  component 
of  geomagnetic  noise,  work  has  bean  done  on  estimating  the 
future  activity  of  aicropulsations,  notably  by  Fraser-Smith 
in  the  case  of  Pci  pulsations  [Ref.  42,  43].  By  extending 
the  predict  ion  technique  for  Pci  pulsations  to  the  Pc2,  Pc3, 
and  Pil  pulsations,  the  occurrence  of  geomagnetic  micropo"- 
sation  activity  in  the  MAD  band  might  be  predicted. 
Combining  this  prediction  with  real-time  solar  flare  infor¬ 
mation  should  give  the  capability  to  disseminate  real-time 
and  estimated  future  HAD  index  values  to  fleet  users. 


Previous  studies  have  led  to  the  conclusion  that  the 
currently  used  geomagnetic  indices  ace  not  accurate  measures 
of  geomagnetic  noise  in  the  HAD  band  pass.  Experimental 
equipment  has  been  utilized,  and  computer  software  written 
in  order  to  confirm  this  conclusion,  and  to  develop  a 
replacement  means  of  evaluating  HAD  geomagnetic  noise. 

A.  EQUIPMENT  CONFIGURATION 

Experimental  equipment,  acquired  as  part  of  the  Naval 
Postgraduate  school  geomagnetics  research  program,  has  been 
utilized  in  the  effort  to  develop  a  usable  MAD  index.  This 
equipment  is  in  use  in  other  projects  of  the  geomagnetics 
research  group.  The  sensors  and  associated  equipment  are  set 
up  for  remote  site  operation  with  system  monitoring  and  data 
analysis  located  at  the  Naval  Postrgraduat e  School. 
Descriptions  of  the  data  collection  system  and  data  analysis 
system  follow  below* 

1 .  £a.la  qoljegtion  System 

The  data  acquisition  system  illustrated  in  Figure 
5.1  reveals  r.he  following  major  components: 

-coil  antenna  sensors  (3) 

-preamplifiers  (3) 

-signal  conditioners  (amplifiers)  (3) 

-pulse  code  modulation  system  (1) 

-radio  transmitter  (1) 

-radio  receiver  (1) 

-instrumentation  tape  recorder  (1) 
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PREAMP 


Figure  5.1  Data  Collection  System 


PREAMP 


a.  Coil  Antenna  sansors 

Three  coil  sensors  are  used  in  this  system. 
Each  sensor  is  a  self-supporting,  continuously  wound,  non- 
center-tapped  coil  antenna  manufactured  by  Sima  Engineering, 
Palo  Alto,  California,  from  about  5460  turns  of  18  gauge 
copper  magnet  wire.  The  coils  weigh  approximately  50  kg  each 
with  dimensions  as  depictad  in  Figure  5.2.  The  dimensions 
of  the  sensor  are  constrained  by  the  dimensions  of  the 
largest  glass  sphere  that  is  commercially  available.  These 
spheres  are  used  to  enclose  the  coils  during  underwater 
experiments.  The  coil  resistance  is  120  ohms  and  its  self¬ 
inductance  is  approximately  9. 31  henries.  The  three  coils 
are  mounted  orthogonaly  on  a  nonmagnetic  frame  (Figure  5.3). 


Figure  5.2  Sensor  Dimensions 


b.  Preamplifier 

The  preamplifier  used  was  the  model  13-10A  low 
noise  ELF  amplifier  manufactured  by  Dr.  Alan  Phillips  of  SRI 


Figure  5.3  Sensor  Mounting  Block 


International.  The  final  stage  of  the  amplifier  contair.3  an 
active  low-pass  filter  which  provides  a  sharp  cutoff  for 
frequencies  above  20  Hz.  The  overall  preamplifier  gain  for 
inputs  of  less  than  2.5  millivolts  is  60  dB. 


f  L 


c.  signal  conditio  ners 


The  signal  conditioners  receive  the  analog 
signals  frost  the  coil  prea  aplifiers,  amplify  them  on  the 
order  of  30  dB,  and  limits  signals  with  peak  amplitudes  of 
7.5  volts  from  entering  the  pulse  cole  modulation  system. 

d.  Pulse  Code  Modulation  (PCM)  system 

The  pulse  coda  modulation  system  chosen  for  use 
is  one  designed  and  manufactured  by  Or.  Bobert  Lowe, 
Lowecom,  Inc.  The  system  features  15  channel  analog  input 
capability  and  offers  selectable  sampling  rages  of  2,  4,  8, 

Id,  32,  64,  and  128  samples  per  second.  By  appropriately 

jumpering  the  analog  input  pins,  the  sampling  rate  may  be 
increased  by  a  factor  of  5. 

The  PCM  system  incorporates  a  crystal  oscil¬ 
lator,  and  associated  CMOS  integrated  circuitry  to  develop 
the  clocking  pulses,  and  a  16  channel  CMOS  analog  multi¬ 
plexer,  a  16  channel,  12  bit  CMOS  analog  to  digital 
converter  and  associated  circuitry  to  provide  the  pulse 
coding.  The  crystal  clock  oscillator  operating  at  a 
frequency  of  24.576  kHz  produces  a  square  wave  output  with  a 
loss  rate  of  1  bit  in  10*.  The  clock  pulses  gate  the  analog 
multiplexer,  analog  to  digital  converter  and  associated 
follow  on  circuitry  that  form  the  pulse  code  words.  The 
basic  output  is  a  Bi-phased  pulse  coded  signal. 

The  data  is  organized  in  frames.  Each  frame  is 
headed  by  a  sync  code  word  which  is  followed  sequentially  by 
the  pulse  coded  samples  from  PCM  channels  1  through  15.  The 
sync  code  word  is  a  pulse  coded  digital  word  with  a  decimal 
value  between  0  and  4095.  This  word  is  preselected  and 
hardwired  on  the  circuit  board.  This  cede  word  is  essential 
to  the  decoding  process. 


In  the  initial  operation  of  this  system,  a 
sampling  rate  of  32  samples  per  second  was  utilized.  Only 
one  sample  per  coil  per  frame  was  analyzed,  making  use  of 
PCM  channels  2,  3,  and  4  only. 

e.  Transmission  and  Becording 

After  the  data  has  been  PCM  encoded,  it  is 
transmitted  by  a  VHP  radio  link  back  to  a  receiver  located 
at  the  Naval  Postgraduate  School,  where  it  is  currently 
recorded  on  an  instrumentation  tape  recorder  for  later 
analysis. 

2.  ha&iisis  saai£JB 

Currently,  the  recorded  PCM  data  is  played  back  into 
a  PCM  decoder  and  associated  equipment  which  generates  a 
nine-track  800  BPI  computer  tape  containing  the  decoded 
sensor  data.  This  computer  tape  is  th’en  input  into  the 

Naval  Postgraduate  School  IBM  3033  computer  for  analysis.  It 
is  in  the  computer  software  that  the  sensor  system  transfer 
function  is  applied,  spectral  analysis  performed,  and  the 
MAD  index  generated. 

B.  DATA  ANALYSIS  SOFTWARE 

As  was  noted  earlier,  a  mainframe  computer  was  utilized 
to  perform  the  spectral  analysis,  apply  the  transfer  func¬ 
tion,  convert  to  power  spectral  density,  plot  the  PSD  and 
generate  the  SMS  noise  amplitude  MAD  index  (using  equation 
4-1).  This  program  is  written  in  FORTRAN  IV  and  is 
discussed  in  brief  below.  A  copy  of  the  program  can  be 
found  in  Appendix  A. 

The  main  program  is  divided  into  sections  which  perform 
the  following  functions: 

-Data  input 

-Fourier  analysis  of  time  series  data 


-Application  of  system  transfer  function 

-Projection  of  field  components  onto  total  field  vector 

-Data  averaging 

-curve  fitting  and  calculation  of  NAD  index 
-Calculation  and  plotting  of  poser  spectral  density 

1 •  Data  Input 

Data  input  is  accomplished  with  the  aid  of  a  subrou¬ 
tine  package  supplied  by  Dr.  Tim  Stanton  of  the  Naval 
Postgraduate  school  Department  of  Oceanography.  Bis  subrou¬ 
tine  (called  ' subroutine  RDM  serves  as  a  fortran  'READ' 
statement,  talcing  the  PCN  data  off  the  computer  tape  and 
converting  it  into  integer  format  with  a  value  between  0  and 
4095.  The  data  input  section  of  the  main  program  takes  this 
integer  value  and  converts  it  to  a  •REAL’  number  and  normal¬ 
izes  it  to  represent  a  voltage  value  between  -5.0  and  +5.0 
volts.  This  section  also  sorts  the  input  data  matching  the 
PCN  channel  to  the  data  array  representing  the  appropriate 
coil. 

2.  zaasisE  A&iizsia 

The  time  series  data  is  next  converted  to  the 
frequency  domain  by  utilizing  a  subroutine  (called  '  FOORT') 
which  performs  a  Past  Fourier  Transform  (FFT) .  The  subrou¬ 
tine  is  one  available  to  users  at  the  Naval  Postgraduate 
School  and  utilizes  the  Cooley-Tukey  FiT  algorithm.  Further 
information  about  this  subroutine  can  be  found  in  the 
program  listing  in  Appendix  A. 

3.  JtefilisaiAaa  si  lii-asAas  Eua&i&ii  ini  Isi&l 
Eiaissiiaa 

The  nazt  section  of  code  applies  the  system  transfer 
functions  for  the  three  coils  to  the  frequency  domain  data. 
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The  transfer  functions  are  given  as  straight  line  segments 
which  were  found  by  least- squares  approximation.  The  data 
enters  this  section  in  amplitude  units  of  volts,  and  is 
converted  into  nanoTeslas  by  the  transfer  function. 

Following  the  application  of  the  system  transfer 
function  to  the  coil  data,  the  program  next  calculates  the 
projection  of  this  data  onto  the  earth’s  magnetic  field 
vector  (total  field  projection).  This  is  done  by  first 

applying  the  local  magnetic  variation  (declination)  to  the 
north-south  (X)  and  East-West  (T)  coil  information  to  deter¬ 
mine  the  horizontal  field  component.  The  local  magnetic  dip 
(inclination)  angle  is  then  used  to  project  the  vertical 
(Z-coil)  and  horizontal  field  fluctuations  onto  the  total 
field  fluctuation. 

4.  Eata  Aisiaaias 

The  previous  program,  sections  exist  inside  of  a 
do- loop  which  enables  the  analysis  of  a  long  period  of  data 
without  a  prohibitive  need  for  storage  apace.  This  loop 
includes  accumulator  arrays  for  each  field  component  and  the 
total  field  projection.  The  fluctuation  data  is  converted 
into  power  prior  to  storage.  This  is  done  by  talcing  the 
magnitude  of  the  fluctuation  component,  dividing  by  the 
number  of  sample  points,  and  then  squaring  the  value.  After 
the  program  passes  through  the  averaging  loop  for  the  last- 
time,  the  arithmetic  average  is  taken  for  each  frequency 
point  on  the  arrays.  At  his  stage  the  power  spectrum  is 
multiplied  by  the  sample  period  to  determine  power  spectral 
density. 

5.  MS  Sal  calami,  on, 

The  next  section  computes  the  RMS  MAD  noise  index 
previously  discussed.  A  polynomial  curve  fit  is  performed  on 
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tha  total  fiald  PSD  using  an  available  subroutine  (*  CHBFT  ’ ) . 
The  resulting  polynomial  is  then  integrated  over  the  limits 
of  the  various  1S3-81  band  pass  settings. 

fiiaiiiac  £a«ai  izis&nl 

Plots  of  the  power  spectral  density  of  each  of  the 
field  components  and  the  total  field  pro j action  are  gener¬ 
ated  in  the  last  section  of  the  program.  This  is  done  by 
converting  the  fluctuation  power  spectral  density  to  deci¬ 
bels  (dB)  referenced  to  1  nanorestia3  per  Hertz.  A  Versatec 
plotting  subroutine  (•PI.OTP1)  is  then  called  to  actually 
generate  the  plots. 


C.  INITIAL  SYSTEM  OPEBATIO  N 

The  NPS  MAD  index  systam  was  initially  placed  into  oper¬ 
ation  with  the  coil  sensors  located  in  the  La  Mesa  village 
housing  area  near  the  Naval  Postgraduate  School,  Monterey, 
California.  System  checkout  was  accomplished  in  June,  1982. 
The  full  system  was  placed  into  operation  on  25  July  19S2 
and  18  August  1982  in  conjunction  with  similar  measurements 
taken  on  the  floor  of  Monterey  Bay. 

The  MAD  index  output  and  power  spectral  density  plots  of 
the  total  field  fluctuation  for  25  July  1982,  1237-1406 
local  (2037-22062),  and  18  August  1982,  0121-0250  local 
(09 2 1 -10 50 Z)  and  0507-0636  local  ( 1307-14362)  are  shown  in 
Tables  IX,  X,  XI,  and  Figures  5.4,  5.5,  and  5.6 
respectively. 
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MAD  Noise  Q921-1050Z ,  18  AUG  82,  Monterey,  CA 
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Figure  5.5  PSD  18  AOG  82,  0921-1050Z,  La  Mesa  Village 
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TABLE  XI 
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MAD  Noise  1307-1436Z ,  18  AUG  82,  Monterey,  CA 
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VI.  CONCLUSION  £££  RECOMMENDATIONS 


A.  CONCLUSION 

Analysis  of  available  information  regarding  current 
geomagnetic  indices  and  the  actual  level  of  geomagnetic 
noise  in  the  MAD  bandpass  indicates  that  the  currently  used 
indices,  the  K  and  Ap  indices,  are  not  valid  for  MAD  opera¬ 
tions.  It  is  therefore  desirable  to  derive  a  new  index 
which  more  accurately  represents  the  geomagnetic  noise  at 
frequencies  of  interest  in  MAD  operations. 

B.  RECOMMENDATIONS 

A  tentative  index  was  developed,  tested,  and  sample  data 
obtained.  The  data  analysis  for  the  preliminary  system  was 
accomplished  using  a  mainframe  computer.  It  is  recommended 
that  further  work  include  improving  the  index  and  the 
setting  up  of  a  system  to  on-line  decode  the  incoming  data, 
and  utilizing  a  desk  top  computer  such  as  the  HP984S  to 
enable  real  time  determination  of  the  MAD  noise  index. 

The  spatial  coherence  of  MAD  geomagnetic  noise  should 
also  be  investigated,  with  possible  application  to  noise 
cancellation.  Additionally,  the  feasibility  of  the  prior 
estimation  of  geomagnetic  noise  should  be  evaluated, 
possibly  by  using  the  prediction  methods  proposed  by 
Praser-Smita  [Ref.  42,  43]. 
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APPENDIX  A 

MAD  INDEX  DATA  ANALYSIS  SOFTWARE 
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